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Abstract 

We show that many well-known signal transforms allow highly ef- 
ficient realizations on a quantum computer. We explain some elemen- 
tary quantum circuits and review the construction of the Quantum 
Fourier Transform. We derive quantum circuits for the Discrete Co- 
sine and Sine Transforms, and for the Discrete Hartley transform. We 
show that at most 0(log^ A^) elementary quantum gates are necessary 
to implement any of those transforms for input sequences of length A^. 



§1 Introduction 

Quantum computers have the potential to solve certain problems at much 
higher speed than any classical computer. Some evidence for this statement 
is given by Shor's algorithm to factor integers in polynomial time on a quan- 
tum computer. A crucial part of Shor's algorithm depends on the discrete 
Fourier transform. The time complexity of the quantum Fourier transform is 
polylogarithmic in the length of the input signal. It is natural to ask whether 
other signal transforms allow for similar speed-ups. 

We briefly recall some properties of quantum circuits and construct the 
quantum Fourier transform. The main part of this paper is concerned with 
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the construction of quantum circuits for the discrete Cosine transforms, for 
the discrete Sine transforms, and for the discrete Hartley transform. 



§2 Elementary Quantum Circuits 

The quantum computation will be done in the state space of n two-level 
quantum systems, which is given by a 2"-dimensional complex vector space. 
The basis vectors are denoted by where x is a binary string of length n. 
The basic unit of quantum information processing is a quantum bit or shortly 
qubit, which represents the state of a two-level quantum system. 

A quantum gate on n qubits is an element in the group of unitary matrices 
W(2"). There are two types of gates that are considered elementary: the XOR 
gates (also known as controlled NOTs) and the single qubit operations. 

The controlled NOT gate operates on two qubits. It negates the target 
qubit if and only if the control qubit is 1. Suppose that x = 6„ . . . 61 is a 
string of n bits, then 



CNOT,,i \x) 



\y) if fcc = l 
\x) if 6c = 



where y is the bitstring obtained from x by negating the bit ht. 

A single qubit gate acts on a target qubit at position t by a local unitary 
transformation 

lan-t (g) ?7(g> l2t-i, U eU{2). 

It will be convenient to describe the quantum circuits with a graphical 
notation put forward by Feynman. The circuits are read from left to right like 
a musical score. The qubits are represented by lines, with the most significant 
bit at the top. Figure |l| shows the graphical notation of the elementary gates. 
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Figure 1: The Feynman notation for the single qubit gate (I2 ® U) and for a 
controlled NOT operation I62 ^'1) = [62 62 © 61). 

A multiply controlled NOT is defined as follows. Let C be a subset of 
[l..n] not containing the target t. Then 



CNOTc.t \x) 



\y) if bc = I for all c G C 
\x) otherwise 



where \y) is defined as above. Several controlled NOT operations in a se- 
quence allow us to implement the operation \x) = |x + 1 mod 2"). Note 
that 0{n) elementary gates are sufficient to realize a multiply controlled NOT 
operation on n qubits, assuming that an additional scratch qubit is available. 
Therefore, at most O(n^) elementary gates are necessary to implement the 
shift operation P„. 
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Figure 2: Shift 

The state of two qubits can be exchanged with the help of three controlled 
NOT operations: 

SWAPfc,/, = CNOT;,,fcCNOTfc,;,CNOT;,,fc. 

It follows that any permutation of the n quantum wires can be realized with 
at most 0{n) elementary quantum gates. 

A more detailed discussion of properties of quantum gates can be found 
in [H . We will discuss the construction of the discrete Fourier transform in the 
next section. In particular, we will show the classical dataflow diagram and 
the corresponding quantum gates to further illustrate the graphical notation. 



•3 Quantum Fourier Transform 



The discrete Fourier transform of length = 2" is defined by 

1 



jk 



J j,fc=0,...,Ar-l 



where uj = exp{2Tci / N) with = —1. Recall the recursion step used in the 
Cooley-Tukey decomposition: 



-N/2 



(Fa ® 1^/2) 



where T; 



N/2 



diag(l, a;, uj"^ 



^/"^ ^) denotes the matrix of twiddle fac- 



tors, and Tjv denotes the permutation given by Tjv \xh) = \hx) with x an 
n — 1-bit integer, and h a single bit. 

We note that the implementation of F2 is a local operation on a single 
quantum bit. The recursion suggest four different parts of the implemen- 
tation of Fourier transforms of larger length. The matrix {F2 ® '^N/2) is a 
single Hadamard operation on the most significant qubit. We would like to 
emphasize that this single quantum operation corresponds to a full butterfly 
diagram. 

The implementation of the twiddle matrix is more complex. Notice 
that T/v/2 can be written as a tensor product of diagonal matrices Lj = 
diag(l,cc;^^ ^) in the form 



T, 



N/2 



Ln-l ® . . .® L2® Li. 



Thus, 1 N/2 ® T]si/2 can be realized by controlled phase shift operations. Fig- 
ure ^ shows the implementation of the two operations discussed so far. 
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Figure 3: For length N=8, only three qubits are necessary. The circuit on 
the left implements {F2 ® '^N/2) and the other realizes the twiddle matrix I4 © 
diag(l,a;,a;^, w^). 



It remains to discuss the other two operations in ([^). The operation 
(I2 ® Fn/2) means that an implementation of the discrete Fourier transform 
of length N/2 is used on the least significant (ra — 1) bits. The operation T^r 
is a permutation of quantum wires. We can combine all the permutations 

TAr(l2 ® Tn/2) ■ ■ ■ (l7V-^2 ® T4) 

into a single permutation of quantum wires. The resulting permutation is 
the bit reversal, see Figure ^ The classical and quantum implementation 
of the discrete Fourier transform of length 8 are compared in Figure 5. We 
observe that the butterfly diagrams find simple realizations but the twiddle 
matrices require more elementary quantum gate operations. 
The complexity of the quantum implementation can be estimated as follows. 
If we denote by R{N) the number of gates necessary to implement the DFT 




Figure 4: The bit reversal permutation is given by I2 T4 followed by Tg. 

of length N = 2" on a quantum computer, then equation (|l]) implies the 
recurrence relation 

R{N) = R{N/2) + 0{\ogN) 

which leads to the estimate R{N) = 0(log^ N). 

Shor's factoring algorithm relies on the quantum Fourier transform in a 
fundamental way. For more details on Fourier transforms and their general- 
izations to nonabelian groups, see 0, H- 



§4 Quantum Cosine and Sine Transforms 

We derive quantum circuits for discrete Cosine and Sine transforms in this 
section. The main idea is simple: reuse the circuits for the discrete Fourier 
transform. 

The discrete Cosine and Sine transforms are divided into various fami- 
lies. We follow p and define the following four versions of discrete Cosine 
transforms: 



,^111 



1/2 

1/2 
1/2 
1/2 



IJTl 
Ki COS — — 



ki cos 
ki cos 
ki cos 



ij=0,...,Af 

i{] -t- l/2)7r 



J i,j=f),...,N-l 



N 

[i + l/2)j7r 

^ J i,j=0,...,Af-l 

:« + l/2)(j + l/2)7r" 



N 



i,j=0,...,N-l 



where fcj := 1 for z = 1, . . . , — 1 and k^ := 1/V2. The numbers ki ensure 
that the transforms are orthogonal. The discrete Sine transforms are defined 




Figure 5: Dataflow graph of the DFT of length 8 and the corresponding quantum circuits. 



by 



ol 



N 



Oil 

On 



cm 



QlV 

On 



1/2 



ki sin 



1/2 r 
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where the constants ki are defined as above. Notice that (resp. S'™) is 
the transpose of CjJ (resp. 5"), hence it suffices to derive circuits for the 
type II transforms. 

It is well-known that the trigonometric transforms can be obtained by 
conjugating the discrete Fourier transform F2N by certain sparse matrices. 
We refer the reader to Wickerhauser for more details on the decompos- 
tions. 



DCTi and DSTi. We derive the circuits for the discrete Sine and Cosine 
transforms of type I all at once. Indeed, the DSTj and DCTi can be recovered 
from the DFT by a base change 



■ F2N ■ Tn — C]^ © iS'^, 



(2) 



where 



/I 



N 



V2 



V2 



1 

V2 
1 

V2 



V2 

i 

'V2 



\ V2 V2 / 

It is straightforward to check that @ holds, see Theorem 3.10 in [0]. Since 
we already know efficient quantum circuits for the DFT, it remains to find 
an efficient implementation of the base change matrix Tn- 

It will be convenient to denote the basis vectors of C^"^^ by \bx), where 
6 is a single bit and x is an n-bit number. The two's complement of an n-bit 
unsigned integer x is denoted by x', that is, x' = 2" — x. The action of Tn 



can be described by 



T^IOO) = |00) T;v|Ox) = -i=|Ox) + ^|lx') 

T^llO) = |10) Tn\Ix) = ;^|0x) -^|la;') 

for all integers x in the range 1 < x < 2"'. Ignoring the two's complement in 
T/v, we can define an operator D by 

Dm = |00) D\Ox) = -i=|Ox) + ^|lx) 

D\10) = |10) D\lx) = 73 lOx) - ^ 

for all integers x in the range 1 < a; < 2". This operator is essentially 
block diagonal and easy to implement by a single qubit operation, followed 



by a correction. Indeed, define the matrix B hv B = —= 
Figure ^ gives an implementation of the operator D. 
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Figure 5: Circuits realizing the block matrix D and the permutation vr. 



Define vr to be the permutation given by a two's complement conditioned 
on the most significant bit vr |0a;) = |0x) and 7r|lx) = \lx') for all n-bit 
integers x. It is clear that Tjv = ttD. The circuits for the permutation vr is 
shown in Figure ^. 

Theorem 1 The discrete Cosine transform C]^ and the discrete Sine trans- 
form Sj^ can be realized with at most 0([og^ N) elementary quantum gates; 
the quantum circuit for these transforms is shown in Figure |^. 

Proof. Let = 2". We note that O(log^A^) quantum gates are sufficient 
to realize the DFT of length 2N. The permutation ir can be implemented 
with at most 0(log^ A^) elementary gates. At most O(logA^) quantum gates 
are needed to realize the operator D. This shows that the DCTi and the 
DSTi can be reahzed with at most 0(log^ A^) quantum gates. The preceding 
discussion shows that Figure ^ realizes the DCTj and DSTj. □ 
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Figure 6: Complete quantum circuit for the DCTj 



DCTiv and DSTiv The trigonometric transforms of type IV are derived 
from the DFT by 



IV 
N ■ 



(3) 



Here Rm denotes the matrix 



/ 1 



R 



N 
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ijj 



UJ 



N-1 



-ILO 



lUJ 



\ 



-lU) 



) 



with UJ the primitive 4A^-th root of unity lo = exp(27ri/4A^). Equation (|^) is 
a consequence of Theorem 3.19 in ^ obtained by complex conjugation. 

Theorem 2 The discrete Cosine transform and the discrete Sine trans- 
form Sj^ can be realized with at most 0{\og^ N) elementary quantum gates; 
the quantum circuit for these transforms is shown in Figure 




Figure 7: Complete quantum circuit for DCTjv 



Proof. It remains to show that there exists an efficient quantum circuit for 
the matrix in equation (^). A factorization of Rn can be obtained as 
follows. Denote by x the one's complement of an n-bit integer x. We define 
a permutation matrix tti by tti |0x) = |0a;) and tti \1x) = \lx) for all integers 
X in the range of < x < 2". Denote by Di the diagonal matrix 

Di = diag(l,co', . . . ,uj^~'^,uj^, . . . ^Zj'^^uJ). 

Then R^ can be factored as 



Rn = 7Vi-D,-\^—\^^ . j (^InJ =7ri-Di-{B®lN). 

Note that B ^ In is a. single qubit operation, and tti can be realized by 
controlled not operations. The implementation of the diagonal matrix Di is 
more interesting. Note that the diagonal matrices of increasing (decreasing) 
powers can be written by tensor products 

Ai = diag(l,u;, . . . = L„ ® ■ ■ ■ (g) L2 ® I^i 

A2 = diag(EJ^-^ ...,uJ,l) = Kn^---(»K2®Ki 

where Lj = diag(l,u;^^ ^) and Kj = diag(cJ^^ Therefore, it is possible 

to write Di in the form Di = (C^In)- (Ai ©A2) with C = diag(l,cJ). The 
circuit for the diagonal matrix Di is shown in Figure M. 
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Figure 8: Quantumcircuit for the diagonal matrix Di. 

The complete quantum circuit for the DCTiv is shown in Figure ^ Note 
that the last three single qubit gates C, B^, and M = diag(e''*/^^, e'^*/^^) 
can be combined into a sing le gate MB^C. □ 

DCTji and DSTh. The implementation of the trigonometric transforms 
of type II follows a similar pattern. Both transforms can be recovered from 
the DFT of length 2N after multiplication with certain sparse matrices, 
cf. Theorem 3.13 in H7|: 



where 



and 



( 1 



N 



1 

71 



V 1 



N 



V2 



V2 



V2 



to) 
'V2 



V2 







V2 



-1 



V2 



/ 



(4) 



and u denotes the 4A^-th primitive root of unity u = exp{27Ti/AN), = — 1. 

Theorem 3 The discrete Cosine transform and the discrete Sine trans- 
form 5*^ can be realized with at most 0{log^ N) elementary quantum gates; 
the quantum circuit for these transforms is shown in Figure |^. 



— H 



4 



271 + 1 



—\Kn~Ln 



t-ri 



c 



T T 

I I 
I I 

1. J. 



-Ki-L 



R 



-1 



B*- J 



-0- 



R 



-1 



Figure 9: Complete quantum circuit for DCTn 



Proof. We need to derive efficient quantum circuits for the matrices Vn and 
U]sf in equation (^. The matrix Vat has a fairly simple decomposition in 
terms of quantum circuits. 



Lemma 4 Viv = 7ri(if ® l^r). 



Proof. It is clear that the Hadamard transform on the most significant bit 
H^Ijy is - up to a permutation of rows - equivalent to Vjy. The appropriate 
permutation of rows has been introduced in the previous section, namely 
TTi |0x) = \lx) and tti = for all < x < 2"-. We can conclude that 
Vat = 7ri(if ® Itv) as desired. □ 

The decomposition of Un is more elaborate. Notice that 



Un |00) = 1 00) Un \0x) = \0x) + \lx') 

V 2 V 2 

Un 111) = (-1) |10) Un \ly) = \0{y + l mod 2")) + \ly) 

for all integers x in the range 1 < x < 2" and all integers ?/in0<y<2" — 1. 
Here and 1 denote the n-bit integers and 2" — 1 respectively. 

Define Dq by Dq\10) = i\lO) and Dq\x) = \x) otherwise. We define 
a permutation 712 by 7^2 \0x) = \0x) and 112 \ lx) = |l(a; + 1 mod 2"')) for all 
integers x in < x < 2". 

Lemma 5 Un = D\T ^ DqT^2- 

Proof. Since D\ |0x) = |0x) and D\ |1x) = u^' |1x), we obtain 
D\Tn\^x)= ^|0x) + ^|1x') 

D\Tn \lx) = |0x) + ^ |lx') 

V2 72 

We have DQn2 |0x) = |0x), DQn2 |lx) = |l(x + 1 mod 2")) for all integers x 
in < X < 2"-l, and Don2 |11) = i |10). We note that (x + 1 mod 2")' = x, 
whence combining DiTn with -Do7'"2 shows the result. □ 
Recall that Tn = ttD. It follows that 

Uj, = n2\DoD')n-'Di. 

The implementation of Di has been described in the section on the DCTiv, 
and the implementation of vr (and hence vr~^) is contained in the section on 
the DCTj. The implementation of vr^^ is also straightforward. It remains to 
find an implementation of -Dq-D*. We observe that 

DoD' \00) = |00) DoD'lOx) = -i= |0x) + ^ |lx) 

DoD' \10) = (-2)110) :DoD*|1x) = 75|0x)-^|lx) 

This can be accomplished by the single bit operation (E)1n followed by a 
multiply conditioned gate J = —j= ^ \ ^ ^ . The full circuit is shown in 
Figure 01. The statement about the complexity is clear. □ 



§5 Quantum Hartley Transforms 

The discrete Hartley transform of length G N is the real N x N matrix 
defined by 



A 



N 



i,j=0,...,N-l 



where the function cas : R — > R is defined by cas(x) := cos(a;) + sin(x), see 
P, p| for classical implementations. The property 



A 



N 



1 - i 



1 +Z 



is easily seen from the definition. We derive a quantum circuit implementing 
Ajsf with one auxiliary quantum bit. 
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Figure 10: Circuit realising a quantum Hartley transform 



Lemma 6 The discrete Hartley transform can be factorized in the form 
shown in Figure \T^. Here R is the unitary circulant matrix 

and H denotes the Hadamard transform. 

Proof. Let F/v be the transformation which effects a DFT conditioned to 
the first bit, i. e., written in terms of matrices we have Fjq = Ijv © -Pat. We 
now show that the given circuit computes the linear transformation |0) \x) ^ 
\Q) An\x) for all unit vectors x G C". Proceeding from left to right in the 
circuit given in Figure [1^ we obtain 

|0)|x) ^ 4(|0) + Il))k) 



^ (|0) + |l))F^|x) 



V2 



I — > 



^|0)F^|a:) + — |0)F^|x) 



-^|0) \^-{l-t)F^ + -{l + z)F^]\x) 

-^|0) \x) + ^\1)F^^An\x) 

i=(|0) + |l))A^ |x) 
|0) A^vk) 



as desired. □ 

Theorem 7 The discrete Hartley transform Ajy can be computed on a quan- 
tum computer using 0{log^ N) elementary operations if we allow one addi- 
tional ancilla qubit. 

Proof. Recall that the discrete Fourier transform Fn can be implemented 
0(log^ A^) operations as shown in Section |§^. The statement follows from 



Lemma ^ since all transformations given there require at most O(log^A^) 
elementary operations. □ 



§6 Conclusions 

We have shown that the discrete Cosine transforms, the discrete Sine trans- 
forms, and the discrete Hartley transforms have extremely efficient realiza- 
tions on a quantum computer. All implementations illustrated an important 
design principle: the reusability of highly optimized quantum circuits. Apart 
from a few sparse matrices, we only needed the circuits for the discrete Fourier 
transform for the implementations. A key point is that an improvement of 
a basic circuit, like the DFT, immediately leads to more efficient quantum 
circuits for the DCT, DST, and DHT. 
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